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We demonstrate efficient heralded generation of high purity narrow-bandwidth single photons 
from a transient collective spin excitation in a hot atomic vapour cell. Employing optical homodyne 
tomography, we fully reconstruct the density matrix of the generated photon and observe a Wigner 
function reaching the zero value without correcting for any inefficiencies. The narrow bandwidth of 
the photon produced is accompanied by a high generation rate yielding a high spectral brightness. 
The source is therefore compatible with atomic-based quantum memories as well as other appli- 
cations in light-atom interfacing. This work paves the way to preparing and measuring arbitrary 
superposition states of collective atomic excitations. 



The last decade has brought about significant advances 
in quantum technology of optical states. Combining the 
techniques of spontaneous parametric down-conversion 
(SPDC), linear optical processing, postselection on mea- 
surement results and optical homodyne tomography has 
led to a plethora of experiments in which interesting op- 
tical states have been prepared and measured [1]. Simi- 
lar accomplishments have been demonstrated in the mi- 
crowave [2] and trapped ion [3] regimes. 

A natural next frontier in quantum state technology 
is to extend methods of quantum state engineering to 
collective spin excitations (CSEs) of atomic ensembles. 
Such extension can find applications in quantum memory 
[4], long distance quantum communication [5], quantum 
logic gates [6], and quantum metrology [7]. In addition, 
quantum engineering within the Hilbert space of atomic 
CSEs is of fundamental interest, as it allows one to ex- 
plore the isomorphism with the Hilbert space of a single 
electromagnetic mode [8]. So far, engineering of CSEs 
has been limited to squeezed spin states [9, 10, 15] and 
the single- quantum state [11, 12]. 

The single CSE quantum can be prepared by heralding 
on detection of a photon that has undergone Raman scat- 
tering from an atomic ensemble, according to the idea of 
Duan, Lukin, Cirac and Zoller (DLCZ) [5]. The Hamilto- 
nian governing the Raman scattering event is identical to 
that of SPDC, leading to the production of a two- mode 
squeezed state of the scattered light and the CSE. While 
DLCZ utilizes only the first-order term of the evolution 
under this Hamiltonian, higher-order terms can be used 
in combination with complex measurements on the scat- 
tered optical mode to produce arbitrary quantum CSE 
states akin to Bimbard et al. [13]. 

Once the desired collective state has been produced, it 
needs to be measured. To that end, the readout stage 
of the DLCZ protocol may be used, in which the CSE 
is converted into the optical domian in a manner sim- 
ilar to readout from a quantum optical memory based 
on electromagnetically-induced transparency [14]. Full 
information about the retrieved optical state, and hence 



about the CSE, can then be acquired using optical homo- 
dyne tomography. An alternative technique of perform- 
ing tomography on atomic CSEs involves off-resonant 
Faraday interactions [15]. 

This outlines an approach to synthesis and measure- 
ment of arbitrary quantum states of atomic CSEs. Here 
we present a proof-of-principle experiment to demon- 
strate the validity of this approach. We produce a her- 
alded single photon from a transient CSE in an atomic 
vapor cell. For the first time for a photon from an atomic 
source, we perform homodyne tomography thereupon, 
obtaining unprecedented uncorrected measurement effi- 
ciency of about 50%, leading to a Wigner function which 
reaches a zero value at the origin of the phase space. In 
this way, our experiment completes the toolbox required 
for complete atomic state engineering. 

Aside from this fundamental aspect, our setup can be 
viewed as a highly-efficient, spectrally bright source of 
single photons for experiments on interfacing quantum 
information between light and atoms. In this context, a 
key requirement is that the photon spectrum be compat- 
ible with the atomic transition in both its bandwidth and 
central frequency. A further crucial figure of merit is the 
spectral brightness measured in the number of photons 
per second per unit bandwidth. Thus far, the majority 
of quantum light sources do not perform well in view of 
these requirements. For example, SPDC-based heralded 
sources, while featuring a high photon production rate, 
also exhibit spectral widths on the order of 10 6 MHz. 
Cavity [16] or waveguide [17] enhanced SPDC sources 
have greatly improved the brightness but both are based 
on filtering a portion of a wider bandwidth and comes at 
the expense of increased experimental complexity. Single 
emitters such as atoms in a cavity [18], quantum dots [19], 
and nitrogen- vacancy centres in a diamond [20] have pro- 
vided on-demand photons albeit with low brightness and 
low collection efficiency. Finally, existing DLCZ-based 
atomic sources of heralded photons, implemented in free 
space [11] and in an optical cavity [12], exhibit sufficiently 
high spectral brightness but relatively poor efficiencies. 



2 



Our source of heralded single photons can be seen as 
a hybrid of SPDC and DLCZ. It is based on four- wave 
mixing (4WM) in hot atomic vapour, which has recently 
been shown to generate high-quality twin-beam optical 
squeezing in a large array of modes [21, 22]. We use 
this twin-beam state to prepare heralded single photons. 
Specifically, we utilize a single, strong, continuous laser 
field that simultaneously pumps both Raman transitions 
associated with the read and write stages of the DLCZ 
protocol [Fig. 1(a)]. In this sense, our system is analogous 
to SPDC: a pair of correlated photons is created simul- 
taneously in the phase-matched directions, which define 
the signal and idler channels. However, the physics of 
our setting is quite different: it utilizes the non- 
linearity in atoms near resonance, while SPDC is based 
on off-resonance m a solid. While SPDC consti- 
tutes conversion of a single pump photon into the sig- 
nal and idler photon, in our process the conversion re- 
quires two photons from the pump field. The required 
third-order optical nonlinearity is sufficiently strong only 
in the neighbourhood of an atomic resonance. Accord- 
ingly, our heralded photon source is intrinsically bright, 
narrowband and produces photons in a well-defined spa- 
tiotemporal mode. While the photons are created off 
resonance to 85 Rb, they are produced directly on reso- 
nance with 87 Rb, and may thus be employed in on- and 
off-resonant protocols. The spectral brightness is among 
the highest observed to date [12, 16, 17], and its mea- 
sured efficiency exceeds that of its competitors by about 
one order of magnitude. 
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FIG. 1. The experiment, (a) The atomic level scheme. Op- 
tical pumping prepares most atoms in the F — 3 hyperfine 
ground state. Raman scattering of the pump leads to the 
creation of a collective spin excitation into the F — 2 state, 
followed by immediate recycling back into F = 3. This pro- 
cess results in emission of a correlated pair of the signal and 
idler photons, (b) The setup. A titanium-sapphire laser is 
passed through a 12-mm cell containing 85 Rb vapour. The 
idler photon passes through spatial and spectral filters and is 
detected on a single photon counting module. Upon a detec- 
tion event, balanced homodyne detection of the signal photon 
is performed. To this end, a phase- stable local oscillator is 
produced by phase-locking a diode laser to the pump laser. 

The experimental setup is diagrammed in Fig. lb. 
The biphotons are generated by passing a single 800- 



mW pump beam of e~ 2 radius 550 fim through a 12- 
mm cell containing 85 Rb vapour at about 100° C. The 
laser is tuned near resonance with the Dl transition at 
A = 795 nm. The signal and idler photons are created at 
a frequency difference equal to the hyperfine ground state 
splitting (3.035 GHz) on each side of the pump frequency 
and are orthogonally polarized to the pump. 

A significant challenge is associated with preventing 
the strong pump field from contaminating the biphoton 
measurements. It is addressed in three ways. First, 
polarization filtering attenuates the pump by 5 orders 
of magnitude. Second, spatial filtering is implemented: 
working in a large Fresnel number configuration allows 
us to achieve phase-matching over a range of angles be- 
tween the signal and idler photons. We choose to work at 
an angle of ±4.7 mrad with respect to the pump which is 
a compromise between phase-matching and spatial sep- 
aration from the pump. Third, frequency filtering of 
the idler channel is implemented. To this end, we em- 
ploy a monolithic spherical Fabry-Perot cavity of a 55- 
MHz linewidth constructed from a standard lens with 
high reflectivity coating on each side, tuned by varying 
its temperature [23]. The cavity is operated slightly off 
the 4WM gain peak in order to filter out the spectral 
region where the idler undergoes Raman absorption, re- 
sulting in uncorrelated photon emission (see, for example, 
[24]). Because of intrinsic spatiotemporal mode filtering 
associated with homodyne detection, filtering in the idler 
channel is unnecessary. In experiments employing pho- 
ton counting, filtering can be accomplished by means of 
a monolithic filter cavity, as in the idler channel. 

Idler photon detection events occur at a rate of « 300 
kHz. Upon each such event, the signal channel is mea- 
sured using a balanced homodyne detector (HD) with 
a 100-MHz bandwidth [25]. The local oscillator for the 
HD must be phase stable with respect to the pump and 
must be at the frequency of the generated signal photons. 
To accomplish this, we lock a diode laser to a frequency 
3.035 GHz blue of the pump beam via an optical phase- 
lock loop [26]. The spatial mode of the local oscillator 
is matched to that of the signal photon by injecting an 
auxiliary laser beam into the idler channel, which mixed 
with the pump inside the cell to generate a classical field 
in the signal mode [27]. 

A quadrature measurement of the heralded state is ob- 
tained from the homodyne current by integrating over its 
temporal mode function ip(t): 

/oo 
q e (t)ip(t)dt, (1) 
-OO 

where qo{t) is proportional to the instantaneous homo- 
dyne detector output photocurrent, and is the phase of 
the local oscillator with respect to the quantum state, 
and the time is measured with respect to the trigger 
event. 
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Knowledge of the temporal mode (1) of the heralded 
photon is essential for its efficient reconstruction [25]. 
This mode is determined by the spectral wavefunction 
of the biphoton produced in 4WM as well as the trans- 
mission spectrum of the filter in the idler channel. The 
former is not known a priori. In order to measure the 
magnitude of the temporal profile, we plot the variance 
of the HD output current as a function of the delay with 
respect to the trigger event [Fig. 2(a)]. For points well 
separated from trigger events, we observe a background 
thermal state, as expected when a two-mode squeezed 
state is traced over one of its modes. In a window of 
approximately 10 ns from the trigger event we note a 
marked increase in variance reaching a peak of 2.3 times 
the vacuum due to the contribution of the heralded sin- 
gle photon to the quadrature noise. To acquire phase 
information of the temporal mode, we observe the auto- 
correlation (q(ti)q(t2)) of the signal, corresponding to the 
density matrix of the heralded state in the time domain. 
The autocorrelation function is of round shape, showing 
high purity of the temporal mode of the heralded photon. 
This observation is further confirmed by the Schmidt de- 
composition of the autocorrelation function. The eigen- 
vector corresponding to the primary eigenvalue yields the 
temporal mode function as shown in Fig. 2(b). We be- 
lieve the side lobes of the mode function to be caused by 
the off-resonant position of the filter cavity. 




FIG. 2. Temporal profile of the heralded photon, (a) The 
magnitude |^(t)| 2 may be inferred by observing the variance 
of the homodyne detector photocurrent as a function of delay 
with respect to the trigger event, (b) In order to infer the 
phase information of the temporal profile, we construct a ma- 
trix corresponding to the autocorrelation of the photocurrent 
as a function of trigger delay. The eigenvector correspond- 
ing to the primary eigenvalue of this matrix (inset) gives the 
temporal weight function ip(t). 

We acquire 10 5 traces of the homodyne current within 
a 180 ns window about the trigger event and integrate 
each trace over the temporal profile determined above to 
produce the ensemble of quadrature values. The quadra- 
ture statistics exhibit no phase dependence, as expected 
for a statistical mixture of Fock states. As seen in 
Fig. 3(a), the observed quadrature probability distribu- 
tion is highly non- Gaussian and shows a well-defined dip 
at the origin characteristic of the single-photon state. If 



the observed state were a perfect single photon, the prob- 
ability to record a quadrature value of zero would vanish, 
but this is not the case, mainly due to the admixture of 
vacuum from losses and mode mismatch. 

In order to determine the density matrix of the de- 
tected state, we apply an iterative maximum-likelihood 
algorithm [28, 29] obtaining an uncorrected single pho- 
ton probability of p tl = 0.488 ± 0.002. The two- 
and three-photon contributions equal p22 = 0.069 an 
P33 — 0.019, respectively. The reconstructed Wigner 
function is shown in Fig. 3(b) displaying a negative value 
of —0.0045 ± 0.0025 at the origin of the phase space be- 
cause the fraction of the odd-number states (primarily 
the single-photon state) is higher than that of the even- 
number states. 




FIG. 3. Experimental results, (a) The marginal quadrature 
distribution for the heralded state (red) and the vacuum (solid 
blue), showing a dip near the origin. The raw set of quadra- 
ture samples is shown below, (b) The Wigner function ob- 
tained from the reconstructed density matrix (upper inset). 
The lower inset shows the cross-section of the Wigner function 
reaching negative values at the origin. 

Since the density matrix gives full information of the 
quantum state, we can calculate common figures of merit 
for single photon sources. We find that in the uncorrected 
reconstructed temporal mode the second order correla- 
tion of the field is g( 2 \0) = (tftfaa)/ (nf = 0.51 < 1 
and the Mandel Q-parameter is Q = —0.32 < 0, evi- 
dencing the nonclassical character of the generated state. 
We note that the difference between the ideal value of 
= and our reported figure is primarily due to 
the higher- number terms. With lower gain correspond- 
ing to a temperature of 69° C, our two-photon compo- 
nent is negligible and is bounded above only by the sta- 
tistical uncertainty of our measurement (0.3%) yielding 
g( 2 )(0) < 0.13, albeit at the expense of reducing the 
single-photon efficiency to 21%. From the quadrature 
variance of the triggered signal (AX t 2 rig ) and the thermal 

background (AX^ ck ), we can also calculate the condi- 
tional cross-correlation between the signal and idler chan- 
nels ff i 2) = «Al t 2 rig > - l/2)/«AX b 2 ck > - 1/2) [16]. At 

(2) 

lower temperatures we find — 24.2 > 1 at zero delay 

(2) 

whereas for higher temperatures we measure = 6.0. 
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From the Fourier transform of the temporal mode, 
we observe the spectral bandwidth to be 39 MHz, cor- 
responding to a spectral brightness of 7, 700 photons 
per MHz per second, comparable to the highest values 
achieved to date [17]. The fundamental limit to spectral 
brightness is set by the requirement that the mean pho- 
ton number of the thermal background state be much 
less than one, to avoid significant overlap of photons. 
The maximum count rate must then be much less than 
the photon bandwidth. In the present work, taking into 
account the efficiency of the single photon detector, we 
are operating on the border of this limit. 

All imperfections characteristic of heralded single- 
photon experiments, such as the linear losses, nonzero 
detection efficiency, imperfect mode matching between 
the local oscillator and the signal [30] exist in the present 
setting as well. Additional imperfections are specific to 
atomic vapour cells, such as uncorrelated Raman scatter- 
ing into the signal and idler channels owing to the popu- 
lation exchange of the hyperfine ground states. This de- 
coherence mechanism is partially overcome by strong op- 
tical pumping by the driving field but this has the detri- 
mental effect of broadening the emission line width and 
introducing significant higher-order terms into the two- 
mode squeezed state of the signal and idler. We therefore 
expect that further improvements to the single-photon 
fidelity can be achieved by incorporating a conservative 
quantity of buffer gas and/or polymer wall coating into 
the vapour cell. This modification will also lead to further 
reduction of the photon bandwidth, since the bandwidth 
of the Raman process decreases with increasing pump 
Rabi frequency. 

To summarize, we have developed a new high spectral 
brightness, high purity source of heralded single photons 
which are compatible with atomic quantum memories. 
The observed uncorrected single-photon efficiency of ap- 
proximately 50% is an order of magnitude higher than 
previously reported figures for atomic photons sources. 
Employing optical homodyne tomography, we recon- 
structed the density matrix of the heralded state and 
found the corresponding Wigner function to be negative. 
Because the photon is obtained from a transient coherent 
spin excitation of the ground state, the experiment can 
be viewed as characterization of the quantum state of this 
excitation, paving the way to a new range of experiments 
on state engineering of atomic ensembles. 
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ful discussions. We acknowledge financial support from 
iCore, AITF, NSERC, and CIFAR. 
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